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I. INTRODUCTION 

The trapping and cooling of fermions has become one 
of the central areas of research within the field of ultra- 
cold atomic gases. Such gases offer the exciting prospect 
of examining the properties of interacting Fermi gases 
with unprecedented flexibility. Due to the existence of 
Feshbach resonances the interactions between atoms can 
be varied almost at will by changing an external magnetic 
field, allowing one to study the intriguing problem of a 
two-component Fermi gas with a unitarity-limited inter- 
action. It has been predicted that the gas at low temper- 
atures consists of a Bose-Einstein condensate (BEC) of 
tightly bound molecules on the molecular side of the res- 
onance where < k-pa <C 1, with a being the scattering 
length and fcp the magnitude of the Fermi wave vector, 
and a Bardeen-Cooper-Schrieffer (BCS) superfluid state 
on the other side of the resonance where < —kva <C 1, 
with an interesting crossover regime in between |lj . Sev- 
eral experimental groups have now reported clear exper- 
imental evidence for a BEC of diatomic molecules on the 
molecular side of the resonance from measurements of 
the momentum distribution of the gas . Experiments 
have also probed the a < side of the resonance using a 
magnetic-field sweep to the BEC side together with radio 
frequency (rf) spectroscopy [3j. 

The study of collective modes is an important spec- 
troscopic tool for probing the many-body dynamics of 
atomic gases 0. Three recent papers report measure- 
ments of the collective modes of a two-component Fermi 
gas ( 6 Li) close to the unitarity limit 0,E3,0]- Although 
the collective mode spectrum of a normal gas in the hy- 
drodynamic regime and that of a bulk superfluid gas are 
identical , it was argued in @ that the damping of the 
modes as a function of temperature should allow one to 



distinguish between a normal gas in the hydrodynamic 
regime and a superfluid gas. 

In the following we shall examine the collective mode 
spectrum of a two-component Fermi gas in its normal 
phase. Our starting point is the Boltzmann-Vlasov equa- 
tion for the fermion distribution function /(r, p, t). The 
effects of atom-atom interactions enter both the collision 
integral and the streaming terms of the equation. Sev- 
eral recent theoretical papers have considered the effect 
of collisions on the collective oscillations of a trapped gas 
of fermions [il ITU UlI IT^ . Our approach differs from that 
of previous authors in that we take into account both the 
energy dependence of the scattering cross section in the 
collision integral and the effects of the mean field in the 
streaming terms. Both these effects are expected to be 
important in the region near the resonance. We deter- 
mine the appropriate relaxation rate from an approxi- 
mate solution to the kinetic equation which is known to 
give very accurate results for the viscous relaxation time 
in the limits of high and low temperatures 0] . We then 
use it to extract the frequency and damping of the col- 
lective oscillations in the trapped gas. 

Since our approach takes full account of the collision 
processes and includes the effects of the mean field, it al- 
lows for a direct comparison between the measured and 
calculated values of the frequency and damping of col- 
lective modes. Our results indicate that the transverse 
oscillations measured experimentally 0, 0j Q are never 
truly hydrodynamic in the normal phase, and it identi- 
fies the regions at high and low temperatures in which 
collisionless behavior is to be expected. Our calculations 
suggest that further insight into the properties of the gas 
near the unitarity limit may be obtained by measuring 
the oscillation frequencies at a fixed magnetic field as a 
function of temperature in a regime where the tempera- 
ture is comparable to or higher than the Fermi tempera- 
ture. 

The plan of the paper is as follows. In Sec. II we in- 
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troduce the viscous relaxation time and derive its tem- 
perature dependence for both a uniform and a trapped 
gas from approximate solutions to the kinetic equation 
for the fermion distribution function. The kinetic equa- 
tion also forms the starting point for the calculation of 
collective mode frequencies and their damping. In Sec. 
Ill we obtain these by taking moments of the equation 
and present our results for a trapped gas of 6 Li-atoms 
in an axially symmetric trap. The effects of interaction 
in the streaming terms are taken into account to lead- 
ing order by adding a mean field to the trap potential. 
For realistic scattering lengths and trap parameters the 
frequencies of the collective modes are close to their val- 
ues in the collisionless limit at high and at low temper- 
atures. In the intermediate temperature regime the gas 
approaches hydrodynamic behavior in the sense that the 
viscous relaxation rate becomes comparable to the oscil- 
lation frequency. In Sec. IV we compare the calculated 
frequencies and attenuation of the collective modes to ex- 
periment, and our main findings are summarized in the 
concluding Sec. V. A detailed account of the moment 
method used for solving the Boltzmann-Vlasov equation 
is given in the Appendix. 



II. VISCOUS RELAXATION TIME 

We shall consider a two-component Fermi gas of atoms 
with mass m in its normal phase. The gas may be uni- 
form or trapped in a potential V(r). We assume that 
the dynamics is described by a semi-classical distribution 
function /(r, p, t) which satisfies the Boltzmann equation 



df ■ df . df 



(i) 



where I is the collision integral. The time-development 
of r and p is given by the equations of motion 



r = v 



_P . 

m 



dV_ 
dr 



(2) 



Since this paper concerns dynamics for which the two 
components of the gas move together, we only need to 
introduce one distribution function referring to a definite 
set of internal quantum numbers, for brevity denoted by 
"spin" with the two values a =t,|- The distribution 
function / always refers to a definite spin value, and we 
have f = f\ = fi- Since we consider s-wave scatter- 
ing, the interaction only involves particles with opposite 
spin. In the present section the streaming terms on the 
left hand side of do not contain any effects of the in- 
teraction, but later, in Sec. Ill, we shall add these as a 
mean field in the equations of motion J2J. 

The viscous relaxation time, which plays an important 
part in the following, is defined in terms of the viscosity. 
Let us briefly recall how one determines the viscosity by 
linearizing the kinetic equation in the spatial derivatives 



of the flow velocity u(r). For simplicity we take the di- 
rection of the flow velocity as our a;-axis and assume that 
it varies in the y-direction, u = (u x (y),0, 0). To calculate 
the viscosity we insert a local equilibrium distribution 
/loc given by 



/loc(p) = /°(e-u-p), 



(3) 



where e = p 2 /2m and / denotes the equilibrium Fermi 
function, and linearize in the gradient du x /dy. Under 
stationary conditions then becomes 



du E v df 
dy v x de 



-I[fl 



(4) 



The shear viscosity r\ relates the momentum current den- 
sity U xy , given by 



n^y 2 



d 3 p 



VyPxf, 



(5) 



to the gradient of the flow velocity according to H xy = 
—rjdu x /dy. The factor of two appearing on the right 
hand side of Q arises from summing over the contribu- 
tions of the two components (| and J,). In order to estab- 
lish the concept of a viscous relaxation time let us make 
a relaxation time approximation to the collision integral 
with a relaxation time r„, which so far is an unknown 
quantity, 



I[f] 



f-f° 



(6) 



The collision integral © yields together with J2J and © 
the viscosity 



(7) 



The total particle density, n to t , of the two components is 



"tot 



d 3 p 



(8) 



The angular integration in momentum space yields a fac- 
tor of 1/15, and the ratio rj/ntot may then after partial 
integration be written in terms of energy integrals as 



ntot 5 / °°63/2/0(l - f o )de 



V 



2 

-T, 



(9) 



Analytic expressions for the integrals occurring in @ 
may be obtained when kT is either small or large com- 
pared to the Fermi energy ep = h 2 kp/2m, where the 
magnitude of the Fermi wave vector, fcp, is given by 
A: F = 37r 2 ntot- When the temperature is much larger 
than the Fermi temperature Tf = ep/k, the equilibrium 
distribution function is / ~ exp[(/i — e)/kT], with ix 
being the chemical potential, and © becomes 



V 
ntot 



kTr„ 



(10) 
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At low temperatures (T < Tp) we obtain from JJjJ that 



n 1 

= -mvpT 

n tot 3 



(11) 



where vp = hkp/m is the Fermi velocity. We shall 
use (jHJl (and the limiting forms (|1C)|) and Qlljl) to de- 
fine the viscous relaxation time in terms of the vis- 
cosity calculated by taking the full collision integral into 
account. The viscous relaxation time at low tempera- 
tures (T <C Tf) would thus be defined in terms of the 
calculated low-temperature viscosity ?yi ow according to 



Tr, = 5ri low /n tot mv F . 



In the presence of the full collision integral /[/], which 
is a functional of the distribution function, we linearize 
the Boltzmann equation in terms of a small deviation Sf 
from the equilibrium distribution by writing /(r, p, t) = 
/°(r,p)+5/(r,p,f) with 

Sf(r,p,t) = /°(r,p)[l - /°(r,p)]3(r,p,t). (12) 

The linearized collision integral becomes a functional of 
$ given by 



/[$] 



d 3 Pi 



^|v- Vl | 



(13) 



where da/dfl is the differential cross section and Q is 
the solid angle for the direction of the relative outgoing 
momentum p' r = (p' — p[)/2 with respect to the rela- 
tive incoming momentum p r = (p — Pi)/2 Q. In the 
present paper, we take the cross section to be given by 
the resonant form 



da a 2 

dQ = 1 + (p r /fr) 2 a 2 ' 



(14) 



In the unitarity limit, where \a\ tends toward infinity, 
the calculated viscous relaxation rate approaches a finite 
value that depends on temperature, since the cross sec- 
tion in this case is determined by the typical value of the 
wave number p r /Ti for the relative motion. 

It should be noted that we neglect in the cross section 
(|14fl any effects of the medium, which can be significant 
at very low temperatures [l4|. 



A. Viscosity of a uniform gas 

We shall now calculate the viscosity of a homogeneous 
gas starting from Eqs. and l|13|) . using a variational 
principle commonly employed in transport theory. The 
Boltzmann equation has the form of a linear, inhomo- 
geneous integral equation, X = H<&, where the inhomo- 
geneous term is X, and H is an integral operator with 
eigenvalues greater than or equal to zero. The use of the 
Schwarz inequality (£/,#[/)($,#$) > (U,H$) 2 , where 
(...,...) denotes a suitably defined scalar product and 



U is a trial function, allows one to put a lower bound on 
the viscosity, which is proportional to (<&, X) = ($, i?$). 
In the present case we choose X to be equal to v v p x and 
define the scalar product according to 



(A,B) 



d 3 p 
(2tt?i) : 



■A(p)B(p)f°(l /°), 



(15) 



where A and B are functions of the momentum. Since 
the momentum-space dependence of the left hand side of 
the Boltzmann equation involves v y p x , it is natural to use 
a trial function proportional to VyJ)x £1 first approxi- 
mation. The viscosity can therefore be approximated by 
the lower bound 



rj = 



2 {X,X) 2 
kf (X, HX) ' 



where the integral operator H is defined as 

1 



/°(l-/°) 



![$]. 



(16) 



(17) 



The explicit factor of two in the numerator of (|16f) arises 
from summing the contributions of the two components. 
The corresponding viscous relaxation time is 



(X,X) 
(X,HXy 



(18) 



which is seen to be independent of the normalization of 
X. The final expression (|18|l is thus an approximate 
expression, obtained by a trial function proportional to 
VyPx , but it is known [13( to differ at high and low temper- 
atures by only a few per cent from the viscosity obtained 
from the exact solution to the Boltzmann equation. 

In the classical limit, T > Tp, and for an energy- 
dependent scattering cross section given by Eq. l(l4^l the 
viscosity, when expressed in terms of the viscous relax- 
ation time t v = rj/ntotkT, is (see, e.g. [T^l 



That 



f I "WT, 1 I 

r n 5y 7r \ m J 



1/2 



(19) 



Here a is an effective cross section, which depends on the 
ratio T/T a , where the temperature T a is defined by 



kT a = 



(20) 



In general, we have 
Ana 2 i" 30 



/ dxx 7 e-*'(l + a?T/T a )- 1 . (21) 
Jo 



For T -c T a we obtain from (|21|l the classical result a — 
Ana 2 , while in the opposite limit, T » T a , Eq. J5TJ yields 



A-na 2 T a 4tt h 2 



3T 



3 mkT 



(22) 
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FIG. 1: The viscosity 77 as a function of temperature for 
fc F |a| = 4.5, in units of ^(Tp) = 5(mkT F /iv) 1/2 /32a 2 , the 
classical value of the viscosity for an energy-independent 
scattering cross section, evaluated at the Fermi temperature 
T — Tp. The inset illustrates the low-temperature T~ 2 de- 
pendence of the viscosity. 

which is seen to be independent of the scattering length a 
and, apart from a numerical constant, equal to the square 
of the thermal De Broglie wavelength. 

At low temperatures, T <C Tp, one expects on general 
grounds that 1/t v oc T 2 due to the restrictions on the 
available phase space caused by the occupied states, the 
so-called Pauli blocking. The magnitude of 1 /r v depends 
on the dimensionless quantity 7 = (kpa) 2 = 2Tp/T a . 
The corresponding variational solution to the Landau- 
Boltzmann equation of a Fermi liquid (see Ref. 0], Sec. 
6.2.1) yields 

1 l-T 2 

-=2n—F( 1 ), (23) 

T r] Til a 

where the function F{^f) is given by the integral 

f 1 x 5 1 
FM = 2 / dx - 5 , (24) 

the variable x being equal to the sine of half the angle 
between the two incoming particle momenta in a colli- 
sion. The function Ffy) decreases monotonically from 
its 7 = value F(0) = 16/15 to its asymptotic expres- 
sion F(7) ~ 4/37 for 7 > 1. 

In the unitarity limit (|o| — > 00) the viscous relaxation 
rate ()23|l becomes independent of the magnitude of the 
scattering length, since F^j) in this limit is proportional 
to 1/a 2 . In general, the calculated relaxation rate tends 
toward a well defined value which depends on temper- 
ature, when the scattering length tends toward infinity. 
The value of 1 /r, ( at unitarity vanishes as T 2 at low tem- 
peratures and as T -1 / 2 at high temperatures. 



In Fig. ^ we plot the calculated viscosity as a func- 
tion of temperature for the value fcp|a| = 4.5. The inset 
shows the viscosity multiplied by T 2 in order to illus- 
trate its characteristic low-tempcrature behavior given 
by (|2S1). Since T a = 0.1T F for this value of k F \a\, the vis- 
cosity at high temperatures is proportional to T 3 / 2 . This 
may be seen by combining the high-temperature relation 
rj = ritotkTTrf with (|19f) and (|22|l . which show that t v is 
proportional to T 1 / 2 at high temperatures, resulting in 
T 3 / 2 . In the case of energy-independent scattering 
(T a ^> Tp) the high-temperature viscosity is proportional 
to T 1 ' 2 . 

B. Viscous relaxation rate of a trapped gas 

In order to apply these results to a trapped atomic 
cloud we now include the trap potential in the equi- 
librium Fermi function. Wc consider the harmonic- 
oscillator potential 

V(r) = j(Lj 2 x x 2 +LJ 2 y y 2 + u 2 z z 2 ). (25) 

The average viscous relaxation rate 1 /r is defined by 

I _ fd?r(X,HX) 

r ~ f<Pr{X,X) ■ 1 ' 

Note that the spatial average of l|18|) is here carried 
out for the denominator and numerator separately. As 
demonstrated in the Appendix, this is the quantity that 
enters as an effective relaxation rate when we take mo- 
ments of the kinetic equation in order to determine the 
frequency and attenuation of the collective modes. The 
calculation of the average viscous relaxation rate pro- 
ceeds as in 0| (see e.g. Eq. (38)), where the correspond- 
ing rate was obtained for bosons above the Bose- Einstein 
condensation temperature, the only modification being 
the change of sign in the equilibrium distribution func- 
tion. The resulting 5-dimensional integrations were car- 
ried out numerically, with varying step sizes until con- 
vergence was achieved. 

The results shown in Fig. [21 and all following figures 
were obtained for a total number N of particles given by 
N = 2.8 x 10 5 , which represents a typical value for the 
experiments on 6 Li reported in [EEI0- We use the trap 
frequencies for the cigar-shaped cloud of 0] , i.e. an axial 
frequency oj z — 2ir x 70 Hz and a transverse frequency 
uj± = 2tt x 1550 Hz, giving an anisotropy ratio equal to 
A = oj z /u ± = 0.045. 

The resulting average viscous relaxation rate is shown 
in Figs. |21 and |3 for two different values of the parame- 
ter fcp I a- 1 , one characterizing the regime of weak coupling 
and the other the regime near the unitarity limit, where 
fcp is the magnitude of the Fermi wave vector in the cen- 
ter of the trap. At low temperatures the relaxation rates 
are proportional to T 2 , and they exhibit in both cases 
a pronounced maximum at a temperature somewhat be- 
low Tp. The asymptotic behavior at high temperatures 
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FIG. 2: The average viscous relaxation rate 1/r divided by 
the transverse trap frequency lu± as a function of tempera- 
ture, for /cf |a| = 0.01. The asymptotic temperature depen- 
dencies are indicated by the dashed lines. Note that the sys- 
tem is highly collisionless, since the maximum value of 1/uj±t 
is about 0.00025. 



FIG. 3: The average viscous relaxation rate 1/r divided by 
the transverse trap frequency lu± as a function of tempera- 
ture, for fcp|a| = 5.5 corresponding to the experiment of |(| 
at a magnetic field of 870 G. The asymptotic temperature de- 
pendencies are indicated by the dashed lines. The dotted line 
is the result obtained in the unitarity limit \a\ — » oo. 



differs in the two cases. When fcp|a| is much less than 
unity, the average viscous relaxation rate decreases as 
1/T at high temperatures. This may seem at odds with 
the fact that for a uniform gas 1/t i( according to H19|) 
is proportional to T 1 / 2 in this limit, since a is indepen- 
dent of temperature. However, the average density in 
a trapped gas is not a constant, but decreases at high 
temperatures in proportion to T~ 3 / 2 , resulting in an av- 
erage relaxation rate proportional to T . When k-p\a\ 
is much greater than unity, the temperature-dependent 
cross section causes the relaxation rate to decrease 
even more strongly, in proportion to T~ 2 . In the unitar- 
ity limit, when \a\ approaches infinity, the average viscous 
relaxation rate approaches a limiting value indicated by 
the dotted curve in Fig. |3] This is further illustrated in 
Fig. 0| where we plot the average viscous relaxation rate 
as a function of 1 / A;p | a | for various temperatures. 

We have normalized in Figs. 121121 and 0] the viscous re- 
laxation rate to the transverse trap frequency lo± used in 
the experiments [f| Q ■ The limiting value of the average 
viscous relaxation rate for \a\ — ► oo is seen never to be 
large compared to u>x, which demonstrates that hydro- 
dynamics cannot be applied to the transverse motion of 
the trapped atomic clouds in the normal phase. In the 
next section we determine the frequency and attenuation 
of the collective modes and obtain results in support of 
this general conclusion. 
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FIG. 4: The average viscous relaxation rate 1/r divided by 
the transverse trap frequency ui± as a function of 1//cf |a| for 
four different temperatures. The parameters T/Tf = 0.03 
and T/Ty = 0.1 correspond to the experimental conditions of 
Refs. |f| and 0, respectively. 



III. FREQUENCY AND ATTENUATION OF 
COLLECTIVE MODES 

We proceed now to extract the dispersion relation of 
the low-lying collective modes by solving the linearized 
Boltzmann-Vlasov equation Q with an appropriate trial 
function The equations of motion are modified to 
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take into account the mean-field potential U given by 

U(r)=gn(v) = g[n°(r) + Sn(v)}, (27) 

where g — 4irh 2 a/m is the interaction constant. The 
density n° denotes the equilibrium density for a single 



spin, that is n° 



'! 



nJ? ot /2, and Sn(r) similarly 



denotes the nonequilibrium change in density for a single 
spin. The effective potential is thus the sum of U and the 
harmonic oscillator potential V(r) given by l|25[) . yielding 
the equations of motion 



P 

— ; p 

m 



d{V + U) 
dr 



(28) 



Let us first establish some useful relations between 
equilibrium quantities. We consider the kinetic equation 
in equilibrium, where the right hand side of Q vanishes, 
and insert the equations of motion l|28|l which result in 



E 



Pi 



or equivalently 



df\ 



f i + mg 



dfi 



dpi 



0, (29) 



E < 

i=x,y,z 



E 



mg 



d[f{i-n\ 



drf \ d[f(l-f )} 



dpi 



(30) 



For convenience we have here introduced the variable 
fi = muiiTi in terms of which the potential is V(r) = 
f 2 /2m. If we now multiply (|3U[) by xp x p 2 and integrate 
over both position and momentum variables, we obtain 



(pIp 



2 i> 



~ x Py) + ~ kT 9 



d 3 r (n ) 2 = 0, (31) 



where (. . .) denotes multiplication by / (1 — / ) and in- 
tegration over the whole phase space, 



= / d 3 j 



d 3 p 



/°(l-/°). (32) 



Using /°(1 — f°) = —(mkT/p)df°/dp, we can calculate 
analytically the in tegr als appearing in (|31|l and obtain 
the virial theorem Hal , 



pot 



^int = 0. 



Here 



Skin — 2 

is the kinetic energy, 



E, 



pot 



d 6 r 



d 3 p p 2 
{2irh) 3 2m 

d 3 p 



(2Trhy 



V(v)f 



(33) 



(34) 



(35) 



the potential energy, while 

E int =g fd 3 r(n ) 2 



(36) 



is the interaction energy. 

If instead we multiply the equilibrium equation l|3U|l by 
the combination x 2 yp y and integrate as before, we obtain 
the equality 



(z 2 pI) 



mg(x 2 y 



dn ^ 



= 0, 



(37) 



which will be used in the Appendix to simplify the ma- 
trix that we diagonalize to obtain the frequencies of the 
collective modes. 

To lowest order in the coupling constant g, the lin- 
earized version of Eq. reads 



E < 

=x,y,z 

mg 



9$ 



dSn df° 



<9$ 

dpi 



- n dfi d P , 



dn° 
-mg— 

dn 



(38) 



where n° as before denotes the equilibrium density for 
a single spin while the corresponding nonequilibrium 
change in the density is 



Sn 



d 3 p 
{2t:Ti)' 



■/ (l-/°)$. 



(39) 



We shall in the following consider modes for which the 
drift velocity u has a spatial dependence given by Ui oc Ti . 
The deviation function $ of a fluid moving with velocity 
u is proportional to u • p. Since acting on u • p with the 
left side of £Q) generates terms like x 2 , p 2 , etc., we follow 
[l?! in choosing the trial function as 

$ = e -«-t J- (aif^ + bihpi + ap 2 ). (40) 

i=x : y ,z 

We insert this ansatz into the kinetic equation 138|) and 
calculate moments by multiplying with the product of 
/°(1 — /°) and any of the terms i 2 , y 2 , . . . ,p 2 appearing 
in $, and subsequently integrating over both r and p. 
The result is a homogeneous set of nine coupled equa- 
tions for the nine coefficients a x ,a y , . . . ,c z and the fre- 
quencies of the collective modes emerge as the roots of 
the determinant. The details of the calculation are given 
in the Appendix for the general case when all three trap 
frequencies arc different. 



IV. RESULTS AND COMPARISON WITH 
EXPERIMENT 

In order to make contact with recent experiments 0, 
0, we consider an axially symmetric trap with lo x = 
Xlu±. We introduce the parameter 



Lo y = uj± and u) 



2E, 



pot 



(41) 
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we find 



-0.2;., 



-0.4- ' 



( =4o;i(l-|; 



(43) 
(44) 



u> 



hd± 



5 + 4A 2 ± 7 ± ^ 4A 4 - A 2 (5 T 7 ) + 2(5 ± 7 ) 



67 



(45) 



-0.6- 



-0.8- 



k p a=-0.5 
. k p a=-1 
. k p a=-2 

k p a=-5.5 



with 7 = (25 - 32A 2 + 16A 4 ) 1 / 2 . For A < 1 the latter are 
approximately given by 



0.5 



1.5 



and 



FIG. 5: The temperature dependence of the parameter £ given 
by (1-411) for different values of kFa. At high temperatures |£| 
decreases as T -5//2 . 



which, as we shall see, determines the sign and relative 
magnitude of the frequency shifts. We shall expand our 
results to first order in £, since our mean-field treatment 
of the interaction in the streaming terms of the kinetic 
equation is only valid when |£| is small compared to unity. 
The temperature dependence of £ is shown in Fig. El In 
accord with our first order treatment of the mean field 
we calculate £ by approximating the equilibrium Fermi 



J hd+ 



J hd- 



10 2; £, 

Y^ 1 + P 



5 zV 20 ; 



(46) 



(47) 



The modes labeled + and — are the transverse and ax- 
ial modes, respectively, which are studied in the experi- 
ments [HE0. 

In the collisionlcss limit we obtain for elongated traps 
(A « 1) 



<4+ = 4^1 



, C V=4^(1-|) 



function, which enters E int as well as -E pot , by its value while for spherical traps (A = 1) 



J poti 

in the absence of interaction. At high temperatures one 
finds from l|41|l that |£| oc T~ 5 / 2 , since the interaction 
energy E- lnii in the classical regime is inversely propor- 
tional to the volume of the cloud (i?; n t oc T~ 3 / 2 ), while 
the potential energy is proportional to the temperature. 

The determinant of the matrix, which is derived in the 
Appendix, has the form of a polynomial in the frequency 
uj. The vanishing of the determinant yields the following 
equation 



^ 2 1+ = 4^ 2 (i+i) 



= 4o,i(l -t). 



(48) 



(49) 



(50) 



(51) 



= LO [{u? - U>l d )-iLJT(LJ 2 - UJ 2 cX )] 



(42) 



where r is defined by (|26|l . Note that the average viscous 
relaxation rate 1/r depends on temperature as illustrated 
in Figs. El El and E| 

In general, the solutions to l)42|) have a real and an 
imaginary part, lo = Re(cj) + iLm(u>), which determine 
the frequency and the damping of the collective modes, 
respectively. The (purely real) frequencies appropriate 
to the hydrodynamic limit, lot — > 0, are denoted by 
subscript (hd), while those for the collisionless limit, 
lot — > co, carry the subscript (cl). To first order in £, 



These results, valid to first order in £, are in agreement 
with those of Pedri et at |ll2j when expanded to first or- 
der in £ , but our results differ to second and higher order. 
This is understandable since the form of Eq. 11401) is more 
general than the scaling ansatz used in |12|. which in- 
volves six rather than nine parameters. However, since 
our calculation of the frequency shifts caused by the in- 
teraction cannot be trusted beyond first order in £, our 
results are in essential agreement with those of 0| . Our 
work thus extends that of ^3 m the sense that we de- 
termine £ and r as functions of temperature, thereby 
allowing a direct comparison with experiment. 

In Fig. El we plot the calculated frequency as a func- 
tion of temperature for fecial = 5.5, which corresponds to 
the parameters used in [6j, along with their experimen- 
tal values. Since we assume |£| to be small compared to 
unity, we show the mean-field curve only in the temper- 
ature region where |£| is less than 0.5. There is a clear 
discrepancy between our calculated frequency and the 
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Hydrodynamic, ^=0 



0.5 1 1.5 1Cf 2 10" 1 10° 10 1 
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FIG. 6: The calculated frequency of the transverse (+) mode 
as a function of temperature, with and without the mean- 
field correction for values of |£| less than or equal to 0.5. The 
experimental values from Q are indicated with the estimated 
error bars included. 



FIG. 8: The frequency of the axial (— ) mode as a function of 
temperature. 



— no mean field 

— with mean field 
experiment 




0.5 



T/T^ 



1.5 



FIG. 7: The inverse damping rate of the transverse (+) mode 
as a function of temperature, with and without the mean- 
field correction. The mean-field corrected curve is plotted for 
values of |£| less than or equal to 0.5. The experimental values 
from (| are indicated. 



experimental result at T — 1.1 5Tp and further work is 
needed to understand the origin of this. The correspond- 
ing results for the damping, given by the imaginary part 
of the frequency, are shown in Fig.0 In order to compare 
with the experimental data below 0.5Tf in Figs. 16171 it 
is necessary to improve our treatment of the interaction 
effects in the streaming terms of the kinetic equation. 

In Figs. OH and [5] we show results for the real and imag- 
inary part of the frequency of the axial mode. Since 
lo z <C lo±_ , there is a broad temperature region where the 
system behaves hydrodynamically. The damping shows 




FIG. 9: The inverse damping rate of the axial ( — ) mode as a 
function of temperature. 



a double-peak structure that reflects, as temperature is 
lowered, the transition between the different regimes, 
from collisionless to hydrodynamic and back to collision- 
less behavior ^j- The mean- field corrections in Fig. |H1 
are seen to be much smaller than those of Fig.[fjl in agree- 
ment with Eqs. lO and ll4"Tj> . 



V. SUMMARY AND CONCLUSIONS 

Starting from a kinetic equation for the semiclassical 
distribution function we have calculated the viscous re- 
laxation rate which determines the frequency and attenu- 
ation of collective modes. By including interaction effects 
as a mean field in the streaming terms, we have solved the 
kinetic equation using a moment method which takes the 
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conservation laws into account and provides an accurate 
account of the damping. Deep in the collisionless regime 
the rate of attenuation is proportional to the viscous re- 
laxation rate, which is small compared to the oscillation 
frequency. If hydrodynamics is applicable, the viscous re- 
laxation rate must be much larger than the oscillation fre- 
quency, and the rate of attenuation is then proportional 
to the oscillation frequency squared times the viscous re- 
laxation time. Our treatment in the present work applies 
to both limits as well as to the intermediate regime and 
yields results that allow for a direct comparison between 
experiment and theory. 

We have demonstrated that for a normal Fermi gas 
over most of the temperature range studied experimen- 
tally 0, 0, 0] j hydrodynamic theory does not apply even 
at the unitarity limit for the transverse oscillations. It 
would be interesting to study further the behavior of the 
gas very near the unitarity limit, for both negative and 
positive a, by measuring the oscillation frequencies at a 
fixed magnetic field as a function of temperature, thereby 
testing the predicted shifts in frequency and attenuation 
due to the interaction. 
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APPENDIX: MOMENTS OF THE 
BOLTZMANN-VLASOV EQUATION 

In this appendix we provide details of the calculation of 
the 9x9 matrix, from which the frequencies and damping 
rates can be extracted. 

First we state some useful identities involving the mo- 
mentum variables pi and the rescaled position variables 
fi = mu)iTi (i = x,y,z), 



together with 



and 



(pIpI) = T;m 2 kTE kin 



(x 2 p 2 x ) = -mtkTEy 



pot ■ 



From these it follows that 



pot 



E, 



pot 



(pIpI) E kia E pot -3E int /2 !-£' 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



Similar identities hold for the other components of p and 
f . They are valid both in the presence and absence of 
the mean field in the equilibrium Fermi function / . 



The following identity holds when the mean field is 
neglected in the equilibrium Fermi function, 



mg(x 2 n°) = mg 



= m 2 kTg Id 



d 3 p 



d 3 r xn° 



—mkT 



df 
dx 



(n°) 2 +x^n° 
ox 



and implies that 



mg(x n ) 
(x 2 pl) 



= — kTE int (A.5) 



(A.6) 



In order to obtain the matrix determining the collective 
modes we insert (|4U|) into <|38ll and start taking moments 
of (|38|) with i 2 , resulting in 

-tw[(3oi +a 2 + a 3 )(x 2 y 2 ) + (ci + c 2 + c 3 )(x 2 p 2 x )} + 

+ uj y b 2 A 2 + u; z b 3 A 3 = 0. (A.7) 

Using the identities given above and Eq. I|37|). we find 

A, = (x 2 p 2 x ) (x 2 x 2 ) mg(i 3 ^) = -2(3?£) 

ox 



An 



(i 2 pl) 



,~2~ dn \ 

ay 



0. 



The constant A 3 also vanishes for the same reason as 
A 2 . Collisions do not appear in (|A.7|) since the combi- 
nations f 2 (as well as Pipi) are collision invariants, i.e. 
I[f 2 } = 0. We now divide Eq. (|A.7jl by {x 2 p 2 ). Since 
the integral mg(y 2 xdn° /dx) appearing in 1)3 7J1 only in- 
troduces frequency shifts of second order in £, we use 
here the approximation (x 2 y 2 ) / (x 2 p 2 ) ~ 1. This yields 

— iuj[3ai + a 2 + a 3 + ci + c 2 + c 3 ] — 2u> x bx — 0. (A. 8) 

Next we proceed to take moments with xp x and obtain 

-iub l {x 2 pl)+L0 x [2{x 2 p 2 x ){ai-ci)+c l m 2 kTE int +B] = 0, 

(A.9) 

where B is given by 

B = —mg 



5n (f° +s °£ 

- /°) V dx 



-mg(x 2 n°) ^ a>i - ^-kTE int c % 

i=x,y,z i=x,y,z 



Dividing Eq. (|A.9|> by (x 2 p 2 ), we obtain 

-iubi + oj x [(2 - £)(ai - ci) - £(a 2 + c 2 + a 3 + c 3 )] = 0. 

(A.1Q) 

Finally we take moments with p 2 . Unlike ff and fiPi 
the quantities p 2 are not separately collision invariants. 
By exploiting the fact that the sum p 2 + p 2 + p 2 — p 2 is 
indeed a collision invariant we arrive at the equation 

- C) _1 (ai + a 2 + a 3 ) + (3 + 4i/3wr)ci 
+ (1 - 2i/3uT)(c 2 + c 3 )] + 2u x b x = 0, (A.ll) 
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where r is given by i|26|) . Eqs. (| A. 8|> . (| A. 1 0|> and (| A. 11|> may be obtained from the three given above by simple 
form respectively the first, second and third lines of our permutation of the indices, 
system of coupled equations. The other six equations 
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